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ABSTRACT 


A renewal-reward model is developed to predict the 
optimum amount of time that Coast Guard personnel should 
spend investigating a vessel for illicit substances. The 
optimal investigation time is determined with respect to 
three criteria; maximizing the number of arrests, 
maximizing the quantity of drugs confiscated, and 
MinimiZing the quantity of drugs that escape detection. A 


Simulation study indicates that the optimal investigation 


time is very sensitive to underlying distributional 
assumptions. The basic service system model may have wider 
application, i.e., to combat modelling, where it may be 


desirable to investigate a potential target to estimate its 
value before committing limited resources. An adaption of 
the model may also be of help in allocating resources for 


mineral exploration. 
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I. INTRODUCTION 


One duty of the United States Coast Guard is the 
interdiction of drug smugglers. When Coast Guard personnel 
board a suspect vessel to search for contraband, a decision 
must be made regarding the amount of time allocated for the 
investigation. If a careful, but lengthy, search is 
conducted, there is a high probability of finding 
contraband if it is present. However, a lengthy 
investigation incurs a penalty. While Coast Guard 
personnel are searching one vessel, other shipping, which 
could have contraband onboard, is passing through the 
patrol area without being examined. At the other extreme, 
if the suspect vessel is examined in a cursory manner, the 
Coast Guard Vessel (CGV) will be able to stop and search 
more shipping but there is a higher probability of not 
discovering the contraband, if present, because of the 
short investigation time. The purpose of this thesis is to 
develop a renewal-reward model to determine the optimum 


investigation time. 


II. RENEWAL-REWARD MODEL 


Consider the following scenario. All shipping 
encountered by a patrolling CGV is stopped and searched to 
determine if contraband is onboard. After leaving port, a 
CGV searcher a time Si until a vessel is sighted. A 
detected vessel has a probability Pn of being Bad (having 
contraband onboard) and Pg = 1- Pp of being Good (having 
no contraband onboard). A Bad vessel contains the random 
number J = 4(1,2,3,...) units of illicit substance with 
probability by. The time required to find the first unit 
of the contraband, given that J units are present, is 
modelled here as the minimum of J independent identically 
distributed times having distribution function Fz (t); this 
model is illustrative only, and may be altered in various 
realistic directions. If time T is required to locate an 


incriminating unit, then 
P (TE =3) = (LR ee ae 


so, upon removing the condition, 


1-Fr(t) = P{Tot} = SS [1-Fz2(t)]J by 3; j=1,2,3,... (2233) 
jel 

Adopt the following decision rule: Establish a 

predetermined investigation time lL. When a vessel is 


detected, it is stopped and searched. If no contraband is 
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discovered in the investigation time L, the vessel is 
released and the CGV resumes patrol. If any contraband is 
discovered before the end of the investigation time L, the 
vessel is detained and escorted to port for further search 
and investigation. A time period D is required to escort 
the vessel to base. The patrol cycle ends when the 
escorted vessel arrives at the base. 

The long run average reward per cycle, R, can be 
calculated in the following manner [Ross: pp. 279-294]. 

E[(Rc ] 


R = 2 2) 
E(C] 





where: E[Rc] is the expected reward for a patrol 
cycle. 


E(C] is the expected duration of a patrol 
cycle. 
The reward associated with an apprehension can be defined 
in various ways. Three different reward criteria will be 
examined. For each of the three cases, we will determine 
the optimal investigation time L which maximizes the long 


run average cycle reward. 


/ EXPECTED DURATION OF PATROL CYCLE: E[C] 
We must distinguish between the following two cases for 
the investigation time I: 


Ig - an investigation time that results in the release of 
the vessel. 


In - an investigation time that results in the detention 
of the vessel. 


Where 


L With probability Pg + Pp Fr (L) 
T with probability PprFr(L) = @ 


| 
o 
it i 


= & 


The length of a CGV patrol cycle can be represented as: 


Si + I» + D With probability a(L) 
Si + Ig + C' With probability a(L) 


where C' has the same unconditional distribution as 


if the first investigation results in a release 


C ..f0n 


of the 


vessel being searched, the process re-starts (regenerates). 


The expected length of a patrol cycle can be expressed 


aS: 

E(C] = E(S] + (Blts] GWE{D] )a ce (Leen tC bic 
or 

B(c) = EiSL+ ElDla + E[Ip]a + La 


l1-a 


which can be written as 


E([C] = E[D] + E[{Ip] + Bis + L& 


RQ] /R 


Solving for the term E[{I»n} yields 


tFr (dt) 
E {Ip ] Fr (L) 


Integration by parts shows that 


— — L 
| Fr (t)dt = Fr(L)L + | tF (dt) 
0 0 


(2.3) 


(2.4) 


(2.5) 


(2.6) 


(2.7) 


Hence 
Pence er (Ll) 


0 
Lh (28) 
Fr (L) 


Therefore 


E(S] + E[D]PpFr(L) + Po | ea tencle ace 


E[C} = 
Pp Fr (L) 
(2.9) 
which reduces to 
" Fr (t)d 
E[C] = E[D] + E[s]_ + , CU SUCKS es est (2.10) 
Pp Fr (L) Fr (L) Pp Fr (L) 
B. REWARD CRITERIA 1 
Suppose the CGV is rewarded for making an arrest. 
Since only one arrest 18S made per cycle, 
E[Re] = 1 (2388) 
and from equation 2.2 
ee 
= = ERG (2.12) 


Consider the following example 


eee 0) , p= 0.08 

S(t) Set vewec me Nr led 

Mio Sie = Des ©.8 

bj = (Leoneye 2 ts , B = 0.8; j = 1,2,. 

Pn = 0.2 
From equation 2.1 

= 2 e-ut (1-8) 

Fr(t) = S(e-!* )! 2(4-6) 6). =) ees (223) 

j=l 1 - Be-"t 
and 
_—_ eer ies ut 
Frit) =416 - Er ©) =22 —esseeee (25145 
1 —- Be-"t 
From equation 2.10, the expected duration of the patrol 
cycle is 
i-6 1-Be-ut 

1 1/r. ( / o,)n( 7 

jC SS 
—ea- ul —ear- ul 
2 Po( 4 e-u /1-pe-»+) @ e-u /1-pe-»s) 
Pg L 
CS Ge ee (2.15) 


—er- ul 
Pp en 2 1-Be-t] 


and from equation 2.12, the long run average reward per 


cycle can be calculated. 





0 a = el 


1 Az Po (2.16) 
ne iaber wt 
L-e- ut Le 1-Be-uL 
eral —— |At * Az *+ As A2 Pp | — Ln aad] LE UPS LA1 X\2 
1-Be-uL Bu 1-8 


Cr REWARD CRITERIA 2 
The reward assigned to the CGV is the amount of drugs 


confiscated. 
E(Rc] = E(J|T < L] (2).17) 


The probability that a vessel being searched contains J = 


jJ units of drugs and that the drugs are discovered is 


P{(T<L)*(J=j)} = E (1-F2 (L))! |b: (2.18) 
so 
P{T<L}@Fr(L) = = [1 2 (1-F2 (L))4 |by (2.19) 
jel 
and thus 
E = (1-Fe (L))4 |ba 
P{(J=5) | (TsL)} = ————————__ (2.20) 
Fr (L) 
and 
= j|1 = (1-F2 (L))4 |b; (2-21) 
jel 
E[J|TsL] = 
Fr (L) 


using the parameters from the Example in Model 1, 


= j[1 : (e-ut ys | (1-B) Bs 

jel 

E(J|F<L] = ————______ (2.22) 
(Loer ee) a -pe-ut) 


which can be evaluated as in the previous example as 


a-e[ (7 / a8) | : (O° Jersey) | 


E[J | T<b] = 





(22240 
using equation 2.2, the long run average reward is 
a E[J Tsb] 
E[C] 
or (2.243 
1 exe 
igre) ( oo | 


e-~ UL 


iL 
aa) a | 
Se, beeen 
ee |S 
———') = ant 


d2 » ( (—=— Py (—— 
T=Be-¥t E=Beg™ = 1-6 es UE 


D. REWARD CRITERIA 3 





The reward examined in model three is the difference 


between the quantity of drugs confiscated and the quantity 





of drugs that escape detection. The expected reward ina 


cycle is 
Rc = E{J|Tsu] - E{M] 


where M is defined as the quantity of drugs that escape 


detection in a cycle. 


Jrse + Jp With probability PpFr (L) 
JIre + Jo +M With probability Pg + PpFr (L) 


<4 
l 


Jrs = Drug quantity passing through patrol area while CGV 
searches a Bad ship 


Jn = Drug quantity passing through patrol area while CGV 
escorts a Bad ship to shore 


Jrs = Drug quantity passing through patrol area while CGV 
searches a Good ship 


Jo = Drug quantity onboard a ship classified as Good 
that is actually Bad 

If conditional expectations are taken assuming that the Bad 

vessels pass through the region according to a Poisson 


process with rate )\1 Pp, then 


E(M] = \1 PoE(J] (E(In]+E[(D]) Po Fr (L)+(A1 PoE(JJE(Ig ]+E(M] ) 


(Ore RHE) tE (do 1Pp Fr (L) (25) 


or 


ety allan { Ps (I Fr (L)at-Lir ) +8 (D1 PF (L) +L (es +e0F «s)} 
E([M] = ¥ NE 
Pp Fr (L) 


E[Jo ] Pp Fr (L) 
+( ; ~~ 
Pp Fr (L) 


where E[Jo] is defined as 


= j(1-F2 (L))I bs; 
j=1 
E(J|T>L] = _- 
Fr (L) 


Using the distributions defined in the Example for Model 1 





Pb A1 16 1-Be7¥Lt 1l-e-- LPg 
E(M|) Se ee — tsi Sea cenmmaemerS + 
ie | ae aie 1-6 Az (1-Be-4t) Pp 
1-Be-¥l 
(1-B)e7uL 
+ == ne ae (2.277 


(1-Be-¥l) (l-e-¥L) 


and the long run average reward can be expressed as 


R = Equation [2.21] - Equation [2.27 (2.28) 
Equation [2.15] 
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ITI. RESULTS 


The optimum investigation time was determined for the 
above example for each reward criteria using two methods. 
The first method consisted of writing a Monte Carlo 
Simulation of the patrol cycle using the distributions of 
the example and running the simulation using various values 
of the investigation time lL. Statistics were gathered 
during the simulation allowing the calculation of the long 
run average reward R. The investigation time was varied 
from zero to six hours in increments of two tenths of an 
hour. The optimum investigation time was determined by 
graphing the long run average reward asa function of 
investigation time and finding the value of L which 
maximized R. Each cycle was replicated 20,000 times. A 
detailed discussion of the simulation can be _ found in 
Appendix A. 

THe second method consisted of writing a computer 
program for the three equations representing the long run 
average reward for the three reward criteria for the 
example and solving the equations for various values of the 
investigation time lL. As in the case of the simulation, 
the value of L was varied from zero to six hours in 
increments) of two tenths of an hour. The optimum 


investigation time was again determined by graphing the 


aca 


long run average reward as a function of investigation time 
and finding the value of L which maximized R. 
Figures 1,2, and 3 contain the results of both the 


simulation and numerical solution using the distributions 


of the examples fOr reward criteria Lara and 3 
respectively. The solid line represents the analytical 
solution and the circles are the simulation results. The 


maximum long run average reward using criteria 1 can be 
achieved by using an investigation time between 1.4 and 2.2 
hours. The maximum long run average reward using criteria 
2 and 3 can be achieved using an investigation time between 
1.2 and 1.4 hours. It is interesting to note that, using 
the distributions presented in these examples, and 
investigation time exists that maximizes all three reward 
criteria simultaneously. This occurs at 1.4 hours’ for the 
input distributions. We donot know that this state of 


affairs will persist. 
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IV. SENSITIVITY ANALYSIS 


In order to check the robustness of the model with 
respect to the underlying assumptions, the simulation was 
run exactly as before with the exception that Fz(t) is now 
assumed to have a lognormal distribution instead of an 
exponential distribution. Each cycle was replicated 20,000 
times. Three separate cases are examined. In case A, the 
distribution function Fz(t) has the same mean and variance 
as the exponential distribution used in the previous 
examples. In case B, the distribution function Fz(t) has 
the same mean but twice the variance as the exponential 
distribution used in the previous examples. Finally, in 
case C, the distribution function Fz(t) has the same mean 
but four times the variance as the exponential distribution 
used in the previous examples. 

Figures 4,5, and 6 contain a comparison between the 
exponential base case and the lognormal cases mentioned 
above for reward criteria 1,2, and 3 respectively. It is 
readily apparent that the results for the exponential base 
case and lognormal case A, with the same mean and variance, 
are different. Furthermore, within the lognormal family of 
curves, it can be seen by examining case B and C that the 


results vary significantly as the variance increases. 
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In order to better understand this behavior, it is 
helpful to compare the quantiles of T, the time required to 
locate an incriminating unit given that J units’ are 


present, for the four distributions. 


Quantile Base Case Case A Case B Case C 
0.10 0.242 1.604 0.908 0.431 
0.25 Oye a2 2.502 1.467 1 632 
0250 2.465 3.945 2.856 1.681 
Cero 6.504 6.858 5010 3.909 
0.90 13.488 t2otoS 107.3385 11.245 


It can be seen that some distributions, in particular 
the lognormal, cases a, b, and ec, have a far greater 
likelihood of producing large values of T relative to the 
others, such as the exponential distribution used in the 
examples. Since the time required to find the first unit 


of drugs greatly influences the cycle length, distributions 





generating larger values of T will produce significantly 
different results. To further verify the accuracy of the | 
distributions produced by the simulation and the quantiles 


listed above, Appendix B contains the calculated quantiles 


| 
) 


—— 


for lognormal case A produced using analytical distribution 


theory. 
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V. CONCLUSIONS 


It has been demonstrated that a renewal-reward approach 
to modelling the Coast Guard drug interdiction process is 
feasible and that it is possible to determine an optimal 
investigation time. By considering several different 
reward criteria in the model, it is possible to quantify 
the costs with respect to the rest of the reward criteria 
when one criteria is selected and used to arrive at an 
optimum investigation time. 

Of the two methods presented to obtain the optimum 
investigation time, the simulation approach is’ the most 
flexible. The example distributions and assumptions used 
in this thesis were selected because they allowed an 
analytical solution to be found. This allowed a comparison 
of results between the simulation and numerical solutions, 
thus verifying the simulation. As the model assumptions 
are changed and different probability distributions 
incorporated, the simulation can be easily modified to 
reflect these changes whereas an analytical solution may no 
longer be possible. 

As the model is currently formulated, its usefulness is 
questionable due to the sensitivity to the underlying 
distribution of the time required to find drugs onboard a 


vessel given drugs are present. As demonstrated in the 


21 


sensitivity analysis even within a given family of 
distributions, the model is sensitive to the distribution 
parameters. In order for the model to produce realistic 
results, rigorous data analysis must be conducted to 
properly identify this distribution. 

The basic service system model may have wider 
application, e.g., to combat modelling, where it may be 
desirable to investigate a potential target to estimate its 
value before committing limited resources. An adaption of 
the model may also be of help in allocating resources for 


mineral exploration. 
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APPENDIX A 


The simulation is written in Fortran 77. It consists 
of a main program and seven. subroutines. Unatrorm [0,1] 
random numbers are provided by calling procedure GGUBFS 
located in the IBM IMSL single precision library. All real 
variables are computed using double precision to minimize 
rounding error. Output is directed to three separate 
units. Unit 2 contains the calculated rewards for each 


increment of investigation time using reward criteria l, 2, 


and 3. Unit 3 contains the cycle length, drug quantity 
confiscated, and quantity of drugs missed for each 
increment of investigation time. Unit 4 contains detailed 


information regarding the cycle for each increment of 
investigation time. This information includes the total 
number of ships searched, the number of ships’ that are 
good, bad, bad but declared good, and bad identified as 
bad. Also included in unit 4 output is the time required 
to discover the drugs on a ship declared bad, the quantity 
of drugs onboard ships passing through the area while the 
CGV is searching a ship, the quantity of drugs missed due 
to short investigation time, and the quantity of drugs 
confiscated on a bad ship. 

The main program controls the starting and final 


investigation time, investigation time increment width, and 
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the number of cycles per investigation time to be 
Simulated. Cycle averages are computed and output directed 
to the three units discussed above. 

Subroutine CYCLE simulates one patrol cycle and records 
all the relevant statistics during the cycle. Subroutine 
STIME generates the random search times required to find a 
ship from an exponential distribution. Subroutine DTIME 
generates the random times required to escort a ship back 
to base from an exponential distribution. Subroutine 
DRGOQTY generates the random quantity of drugs on a BAD ship 
from a geometric distribution. Subroutine CLASS determines 
the classification of a vessel based on the deterministic 
value Pp. Subroutine OPCOST determines the quantity of 
drugs missed onboard other ships passing through the area 
while the CGV is investigating the current ship. 
Subroutine RTIME generates the random times required to 
find the first unit of drugs given that Jj units are 
present. 

Four versions of this subroutine are listed; one for 
the exponential distribution and the other three for the 


lognormal case A, B, and C distributions. 
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FILE: DPSIM FORTRAN Al 


THIS PROGRAM SIMULATES PATROL CYCLES FOR VARIOUS LENGTHS OF 
INVESTIGATION TIMES AND CALCULATES THE VALUES FOR THE THREE 
REWARD CRITERIA. FUNCTION GGUBFS FROM THE IMSLSP LIBRARY IS 
USED TO GENERATE UNIFORM (0,1) NUMBERS. 


INTEGER I,K,MAXT,NCYCLE,Z 

REAL¥8 L,C,J,M,INC,EC(500),EJ¢500),EMC 500), ITIMEC500),TOEND, 
& R1(500),R20500),R3(500),TS,GS,BS, BSG, BSB,EIB,OP,MISS,EIG, 
& CATCH,ES,ED,BJ, DSEED 

DSEED = 995317.D0 


ANMAND 


L = 0 
TOEND = MAXT/INC 
DO 200 = 1,TOEND 
L + INC 
C 
J 
M 0 
PRINT®, ‘COMPUTING L = *,L 
DO 100 I = 1,NCYCLE 
CALL CYCLECDSEED,L,C,J,M,1TS,GS,BS,BSG,BSB,EIB,OP,MISS, 
& CATCH,ES,ED, BJ, EIG) 
100 CONTINUE 
Z=NCYCLE 
WRITECG,®) L 
WRITEC4,523) TS/Z,GS/Z,BS/Z,BSG/Z, BSB/Z 
MIRITEC4,524) EIB/Z,BJ/Z,0P/Z,MISS/Z,CATCH/Z 
WRITEC4,525) EIG/Z,ES/Z,ED/Z 
WRITECG, ®) 
525 FORMAT C1X,5F9.4) 
524 FORMATC1X,5F9.4) 
sae FORMATC1X,3F9.4) 


TS =0 
BS =0 
BSG=0 
BSB=0 
GS =0 
EIB=0 
MISS=0 
pacar 


oOOraA 


mM 
om 
i | 


© 
“O 
i 
ooooaeo 


EIG= 
ITIMECK) = L 
ECC K) C/NCYCLE 
EJCK) J/NCYCLE 
EMCK) M/NCYCLE 
RICK) LAGU) 
R2CK) EJCKIJ/ZECCK) 
R3CK) CEJCK)-EMCK) DZECCK) 
WRITEC2,19) ITIMECK),RICK),R20K),R35CK) 
WRITEC3,19) ITIMECK),ECCK),EJCK),EMCK) 
19 FORMAT (1X,F6.2,5F9.3) 
200 CONTINUE 
aoe "PROGRAM COMPLETE’ 
D 


SUBROUTINE CYCLECDX,LSTAR,CSTAR, JSTAR,MSTAR,TS,GS,BS,BSG, BSB, 
EIB,OP,MISS,CATCH,ES,ED, BJ, EIG), 

REAL¥8 LSTAR,CSTAR,JSTAR,MSTAR, LAMDAL,LAMDA2,MU, BETA,S,R,QTY, 
& D,LEAK, TIME,PBAD,TS,GS,BSG,BSB,EIB,OP,MISS,ES,ED,BJ,BS,CATCH, 
& EIG, DX 

INTEGER NC 

CHARACTER®4 TYPE 

LAMDA1 = 1.10 

LAMDA2 = 0.3 
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PICE: BPSIK FORTRAN Al 


20 


ZA 


MU = QO. 
BETA = 


CALL STIMECDX,LAMDAI, 5S) 
ES=ES+S 
CALL CLASSCDX,PBAD, TYPE) 
IFC TYPE .EQ. 'GOOD') THEN 
EIG = EIG+LSTAR 
GS=GStl 
QTY=90 
TIME = LSTAR 
CSTAR = CSTAR + TIME + §S 
CALL OPCOSTCDX,LAMDA1, TIME, BETA,PBAD,LEAK) 
MSTAR = MSTAR + LEAK 


OP=OP+L EAK 
GOTO 20 

ELSE IFC TYPE .EQ. 'BAD') THEN 
BS=BStl 
CALL DRGQTYCDX, BETA, QTY) 
BJ=BJ+QTY 


CALL RTIMECDX,MU, QTY,R) 

IFCR .GT. LSTAR) THEN 
EIG=EIG+LSTAR 
BSG=BSG+t1 
TIME = LSTAR 
CSTAR = CSTAR + TIME + S$ 
CALL OPCOSTCDX,LAMDA1, TIME, BETA,PBAD, LEAK) 
MSTAR = MSTAR + LEAK + QTY 
OP=OP+LEAK 
MISS=MISS+QTY 
GOTO 20 

ELSEIFCR .LE. LSTAR) THEN 
EIB=EIB+tR 
BSB=BSB+l 
CALL DTIMECDX,LAMDA2,D) 
ED=ED+D 
TIME = R + D 
CSTAR = CSTAR + TIME + S 
JSTAR = JSTAR + QTY 
CALL OPCOSTCDX,LAMDAL, TIME, BETA,PBAD,LEAK) 
MSTAR = MSTAR + LEAK 
OP=OP+LEAK 
CATCH=CATCH+QTY 
GOTO 21 

ELSE 
STOP ‘ERROR 1° 

ENDIF 


E 

STOP *ERROR 2° 
ENDIF 
RETURN 
END 


SUBROUTINE STIMECDRSEED, INPUT,OUTPUT ) 
REALX%8 INPUT,OUTPUT,RV,DRSEED 
RV=GGUBFSCDRSEED) 

OUTPUT = LOGCRV)I/C-INPUT) 

RETURN 

END 


SUBROUTINE DTIMECDRSEED, INPUT, OUTPUT) 
REAL¥8 INPUT,QUTPUT,RV, DRSEED 

RV = GGUBFS( DRSEED) 

OUTPUT = LOGCRV)/(C-INPUT) 

RETURN 

END 
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moe DPSIM FORTRAN Al 


SUBROUTINE DRGQTYCDRSEED, INPUT, OUTPUT ) 
REAL®8 DRSEED,INPUT,OUTPUT,RV 

RV = GGUBFSCDRSEED) 

OUTPUT = AINTCLOGCRV)/CLOGCINPUT)))+1 
ae 


SUBROUTINE CLASSCDRSEED, INPUT,OUTPUT ) 
REAL¥8 DRSEED, INPUT,RV 
CHARACTER¥®4 OUTPUT 
RV = GGUBFS(CDRSEED) 
TEGRVecce. INPUT) THEN 
OUTPUT = 'BAD' 


Se 
OUTPUT = 'GOOD' 
ihe 


RETURN 
END 


SUBROUTINE OPCOSTCDRSEED, INPUT1, INPUT2,P1,P2,QUTPUT ) 
REAL*8 DRSEED, INPUT1,INPUT2,P1,P2,O0UTPUT,STOR,NSHIP 
OUTPUT = 0 

JO Vser7 C1-P1D 

NSHIP = INPUT2X*INPUT1 

OUTPUT = (CP2¥*NSHIP)XJO 

RETURN 

END 


SUBROUTINE RTIMECDRSEED, INPUT1L,INPUT2, OUTPUT) 
EXPONENTIAL DISTRIBUTION - BASE MEAN & VARIANCE 
REAL¥®8 DRSEED, INPUT1,INPUT2,OUTPUT,RV 

RV = GGUBFSCDRSEED) 

OUTPUT = LOGCRV)/C-INPUT1LXINPUT2) 

RETURN 

END 


SUBROUTINE RTIME1LCDRSEED, INPUT1, INPUT2, OUTPUT) 
LOGNORMAL DISTRIBUTION - BASE MEAN & VARIANCE 
REAL*¥8 INPUT1,INPUT2,OUTPUT,STDV,AVG,PI,LOW,U1,U2,GG, 


& DRSEED 


INTEGER II,FINISH 

GG = LOG(2.0) 

LOW = 999999 

FINISH = INTCINPUT2) 

STDEV = SQRTC(GG) 

AVG = -CLOGCINPUT1))- 0.5%GG 
PI = 3.141592654 

in 662 II] = 1,FINISH 


= GGUBFS( DRSEED) 
U2 = GGUBFSCDRSEED) 
RV = SQRTC-2*¥LOGCUL))¥COSC2*PIXUZ2) *STDEV+AVG 


RVT = EXPCRV) 
IFCRVT .LT. LOW) THEN 


LOW = RVT 
ENDIF 
CONTINUE 
OUTPUT = LOW 
RETURN 


END 


SUBROUTINE RTIME2CDRSEED, INPUT1,INPUT2,OUTPUT ) 
LOGNROMAL DISTRIBUTION - BASE MEAN & 2*VARIANCE 
REAL*%8 INPUT1,INPUT2,OUTPUT,STDV,AVG,PI,LOW,U1,U2,GG, 


& DRSEED 


INTEGER II,FINISH 
GG = LOGC3.0) 
LOW = 999999 


Ja7 | 


FILE: 


662 


662 


DPSIM FORTRAN Al 


FINISH = INTCINPUT2) 
STDEV = SQRTC(GG) 
AVG = -CLOGCINPUT1))- 0.5%GG 
PI = 3.141592654 
DO 662 II = 1,FINISH 
GGUBFSCDRSEED) 
GGUBFSCDRSEED) 
SQRTC-2*XLOGCUL) ) ¥COSC2¥PIXUZ)*XSTDEV+AVG 
RVT = EXPCRV) 
IFCRVT .LT. LOW) THEN 

W = RVT 


AC 
<N 
tb ok Mb 


CONTINUE 
OUTPUT = LOW 
RETURN 

END 


SUBROUTINE RTIME3(DRSEED, INPUT1, INPUT2, OUTPUT) 
LOGNROMAL DISTRIBUTION - BASE MEAN & 4¥VARIANCE 
REAL¥8 INPUT1, INPUT2,QUTPUT,STDV,AVG,PI,LOW,U1,U2,GG, 


& DRSEED 


INTEGER II, FINISH 
GG = LQG(5.0) 
LOW = 999900 
FINISH = INTCINPUT2) 
STDEV = SQRT(GG) 
AVG = -CLOGCINPUT1))- 0.5XGG 
3.141592654 
62 II = 1,FINISH 
GGUBFSCDRSEED) 
GGUBFSCDRSEED) 
SQRTC-2*¥LOGCUL))¥COSC2¥PI RUZ) XSTDEV+AVG 
RVT = EXPCRV) 
IFCRVT .LT. LOW) THEN 
LOW = RVT 


ENDIF 
CONTINUE 
OUTPUT = LOW 
RETURN 

END 


oo 
Oo 
th ok bb ON Oh 
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FILE: TSTNUM FORTRAN Al 


C 
Cc 


100 


THIS PROGRAM COMPUTES THE VALUES OF THE THREE REWARD CRITERIA 
FOR VARIOUS VALUES OF L FOR COMPARISON WITH SIMULATION RESULTS 


REAL PB,LAMDA1,LAMDA2, BETA,MU,INC,ES,ED,FT,FTBAR,A,ABAR, 
& L,TOEND,EC(100),EJ(100),EM(100),R10100),R2¢100),R3¢0100), 
& ITIMEC100),PG,B,BB,E,C1,C2,C3,C4,M1,M2,TEMP,J1,RATE 

eres MAXT,K 


MAXT = 6 
INC = 0.20 
PB = .20 
PG = 1-PB 


LAMDA1L = 1.10 
LAMDA2 = 0.3 
ES=1/LAMDA1 
ED=1/LAMDA2 
BETA = 0.8 

MU = 0.08 


RATE = CLAMDAI1X®PB)/BB 

Meant = MAXT/INC 

DO 100 K = 1,TOEND 
L=L+INC 
E=EXPC-MUXL ) 
TEMP=1-(BXE) 
FTBAR = (BBXE)/TEMP 
FT=1-FTBAR 
ABAR=PBXFT 
A=1-ABAR 
Cl= ES/ABAR 
C2=ED 
C3=(PB/ABAR) XC CBB/ CBXMU) )XLOGCTEMP/BB)-CL¥FT BAR) ) 
CG=(L¥A)/ABAR 
JI=CBB/FT)¥CC1/CBB¥BB)I-CE/CTEMPXT EMP) )) 
ML1=C(C3+ED) XRATE+C CL XAXRATE)/ABAR) 
M2=BBXCE/CTEMPXTEMP) )®¥C1/ABAR) PB 
ITIMECK) = L 


ECCK) = C1+C2+C3+C4 

EJC(K) = Jl 

EMCK) = M1+M2 

RICK) = 17ECCK) 

R2CK) = EJCKIZECCK) 

R3(K) = CEJCK)-EMCKII/ZECCK) 


WRITEC2,132) ITIMECK),RICK),R2CK),R3CK) 
WRITEC3,132) ITIMECK),ECCK),EJCK), EMCK) 
WRITECG,®) L 
WRITEC4$,133) C1,C2,03,C4 
WRITEC4,134) Jl 
WRITEC4,135) M1,M2 
FORMAT (1X,F6.2,35F9.3) 
FORMAT (1X,4F12.5) 
FORMAT (1X,F12.5) 
FORMAT (1X,2F12.5) 
WRITECS, ®) 
CONTINUE 
STOP "PROGRAM EXECUTION COMPLETE’ 
END 
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APPENDIX B 


Recall. that the probability distribution for the 


quantity of drugs on a vessel given that it is bad is: 


by = (1.-— 8)637!3— © — 0263) 12 


Let 
P{T>t} = Fr(t) = SE (1-Fz (t)) Ib; 
j=l 
so 
=a _ (1-8) Fz (t) 
Bx (ty i -pFr (th 
and 
_ , _ (1-B)F2(t) _ _1-F2 (t) 
PME 1-BFz (t) 1-BF2 (t) 
now let 
iSE(te) 2 
D165) 0) (tp ) 


which can be written as 


Fz (tp ) = re 


and 


1 - Fz (tp ) = ue 


Fz (tp) 1-Bp 
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Assuming Fz(t) is lognormal with mean = 2.1792 and standard 


deviation = 0.8326 as in distribution Case A 


Q (is tp - ») jo ed OB, 
fo) 1-Bp 


therefore 


tp = EXP W+ OZp(1-f ) 
“1-fhp 


where Zpi(i-p) 18S taken from tables of the Standard Normal 
1-BPp 


distribution. 


Using the equation for tp, the following table of quantiles 


was generated: 


p tp 
0210 1.65 
OF.25 2.46 
0.50 3295 
0.75 6.78 
0.90 11.99 


It can be seen that these quantile values agree closely 


with those presented in the sensitivity analysis. 
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